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en ' Abstract 

> 

cn . We consider the minimal coupling of two photons to neutral scalar and pseu- 

doscalar fields, as for instance in the case of the Higgs boson and axion, respectively. 
"^ , In this framework, we analyze the photon dispersion relations in the presence of 

vQ \ static and homogeneous external magnetic and electric fields, by taking into account 

^D \ the contribution of the imaginary part of the scalar/pseudoscalar self-energy. We 

,^ ' show that this contribution cannot be neglected when it is of the same order as the 

photon-scalar/pseudoscalar mixing term. In addition to the usual light-like photon 
propagation mode, with a refraction index n > 1, a massive mode with mass of the 
,^ . order of the coupled boson mass can be induced, provided that the external field is 

above a particular critical value. Depending on the values of the external field, pho- 
k> \ ton energy, and mass of the scalar/pseudoscalar particle, the scalar/pseudoscalar 

H ' width could induce a sizeable rate of photon splitting in two photons due to a 

strong resonant phenomenon. This effect has no practical laboratory applications 
for the Higgs physics due to the very large critical external magnetic or electric 
fields involved, for a photon energy of the order of a TeV. However, it can have rel- 
evant consequences in the axion physics or in any other scenario where light neutral 
scalar/pseudoscalar fields have minimal coupling with two photons. 



ix. 



1 Introduction 

Scalar and pseudoscalar neutral particles can have effective gauge invariant coupling with 
the electromagnetic field A^ due to the higher dimensional operators. In particular, in the 
case of a scalar {^ps) and pseudoscalar {^pp) fields, the minimal interaction Lagrangians, 
proportional to the lowest dimensional gauge invariant operators, are given respectively 
by 

^ Fm.F^>5, (1) 
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where F^,y = df^A^, — d^A^ is the electromagnetic field strength tensor, F^j, = \^f_iua(3 F^^ 
its dual, and the mass scales Ks^p parametrize the corresponding couplings. 

Well known examples of such couplings to two photons can be found, for instance, 
in the framework of Higgs boson [1,2] and axion physics [3-5]. The Higgs mechanism, 
responsible for the electroweak symmetry breaking (EWSB) in the Standard Model (SM), 
predicts a heavy scalar particle, the Higgs boson, which up to now has eluded any exper- 
imental test. The Higgs boson mass cannot be predicted in the framework of SM, and it 
ranges from the present experimental lower limit [6, 7] of 114 GeV, up to the theoretical 
upper bound of around 800 GeV [1]. Although the Higgs boson has no tree-level coupling 
with two photons, it acquires such a coupling at 1-loop [2], with the corresponding scale 
As of the order of O(TeV). 

On the contrary, the axion, which arises as a natural solution to the strong CP problem 
[3,4], is a very light pseudoscalar particle which remains undetected due to the very 
weak coupling to ordinary matter. It appears as a pseudo Nambu-Goldstone boson of a 
spontaneously broken Peccei-Quinn symmetry [3] and has a mass itia expected to be of 
the order of itia ~ (9(meV) [8-11]. Moreover, since the corresponding coupling Ap is of 
the order of 10^° — 10^^ GeV if astrophysical constraints are taken into account [11], the 
decay width of the axion into two photons, being proportional to m\/ K^p, is very small, 
implying that the axion is almost a stable particle on cosmic time scales. 

The fact that a scalar field has an effective coupling with two photons, allows the 
'-) -^ ifs conversion mechanism in external magnetic or electric field to work by means of 
the Primakoff effect [12]. A simple way to look at this phenomenon is the following. By 
replacing one of the photon fields in ^psll coupling in Eq.(l) with an external electro- 
magnetic (EM) background field, the resulting (^957 mixing term could be regarded as an 
off-diagonal term of the scalar-photon mass matrix, providing a source for the '-/ ^ ips 
transition. This is not in contradiction with angular momentum conservation, since due to 
the presence of the external field, the SO (3) rotational invariance for the total system ra- 
diation + background is restored. Analogous results hold for the pseudoscalar conversion 
mechanism 7 -^ ifp as well [9, 10]. 



The mechanism of photon-axion conversion in external magnetic field has been in- 
tensively studied in the literature [8-11], especially in the context of astrophysics and 
cosmology [11], as the time evolution of the quantum system of two particle states. Due 
to the very large life-time of the axion, the photon-axion oscillation phenomenon is pos- 
sible on macroscopic distances. However, it is worth stressing that the photon-axion con- 
version, with both fields on-shell, is not possible on strictly homogeneous and constant 
magnetic fields in all space. Indeed, due to the fact that the axion is a massive parti- 
cle, the conservation of the energy and momentum at the axion-photon vertex requires 
inhomogeneity of the external magnetic field in order the reaction to proceed. Therefore, 
the transition amplitude comes out to be proportional to the form factor of the external 
field, inducing a suppression [8-10]. When the scale (L) over which the magnetic field 
is homogeneous, is comparable or smaller than the scale of the inverse of the transfered 
momentum AP ~ m\/{2E), the magnetic form factor could easily be of order 0{1). 
This is indeed the case of optical photons and axion masses of the order of meV, where 
APL ^ 1. On the contrary, when APL ^ 1, a strong suppression is expected from the 
form factor if the magnetic field is homogeneous. However, the depletion induced by the 
form factor can be ameliorated by embedding the system in a plasma, and tuning the 
plasma frequency, which provides an effective mass for the photon, to be equal to the 
axion mass [13]. 

A particular class of experiments [14, 15] have been using the idea of the photon 
regeneration mechanism [9] in order to detect the axion in laboratory. After a Laser 
beam passes through a magnetic field, an axion component can be generated. The original 
incident Laser is then absorbed through a thick wall, but not the axion component. This 
can be then re-converted into photons by passing through a second magnet, the so-called 
phenomenon of light shining from a wall. 

However, indirect effects of the axion coupling with two photons could give rise to the 
modification of photon dispersion relations in the presence of an external constant mag- 
netic field [16]. In particular, these effects would induce the so called 'Faraday rotation', 
i.e. the rotation of the plane of polarization of a plane polarized light passing through a 
magnetic field [16]. The BFRT collaboration [15] has also performed a polarization exper- 
iment along these lines, but no significant deviation has been observed [14]. Recently, the 
PVLAS experiment [17], has reported the first evidence of a rotation of the polarization 
plane of light propagating through a static magnetic field. If confirmed, these results 
might point out the presence of a very light pseudoscalar particle weakly coupled to two 
photons. Alternative experiments have been proposed to check these results [18, 19]. The 
phenomenon of 'Faraday rotation' in a medium has also been analyzed in the context of 
QED [20]. 

In this framework, the aim of the present paper is to investigate the photon dispersion 
relations on static and homogeneous external electric or magnetic fields, by using the 
approach of the effective photon propagator. This last one is obtained after integrating 



out the scalar /pseudoscalar degree of freedom, or equivalently by summing up the relevant 
class of photon self-energy diagrams as done in our paper. Apart from the different 
approach, the new aspect with respect to previous studies mainly consists in the fact that 
we are including the contribution of scalar/pseudoscalar self-energy diagrams, where the 
effects of both real and imaginary parts are taken into account. In particular, we will 
show that the contribution of the imaginary part cannot be neglected. This approach 
would allow to analyze new phenomena not considered before, such as the effects of 
the scalar/pseudoscalar decay width (F) on the absorptive part of the effective photon 
propagator. Nevertheless, one should also expect that the width T corrections could induce 
some effects on the rotation of the polarization plane of the incident photon, although we 
have not considered this aspect in the present work. 

Analogous studies are well known in the framework of QED, where the refraction 
index for photon propagation and absorption coefficient have been determined in the case 
of photon propagating in strong magnetic fields [21-23]. In our case, this task can be 
easily achieved by summing up the full series of leading Feynman diagrams of photon 
self-energy and dispersion relations can then be derived by looking at the poles of the 
effective propagator. While the real part of the scalar/pseudoscalar self-energy (ReS) can 
be re-absorbed in the renormalization of the corresponding mass and wave function, the 
ImS would induce an imaginary part in the effective photon propagator, giving rise to a 
no n- vanishing contribution to the photon absorption coefficient. 

Due to the presence of the external field, the vacuum polarization and the dispersion 
relations will be modified with respect to the case without external sources. This will affect 
only the component of the photon polarization vector interacting with the external field, 
according to the Lagrangian in Eqs.(l), (2), while the other one will remain unaffected 
[16]. In particular, there will be two new propagation modes in addition to the one which 
is not modified. The lightest mode, corresponding to a light-like particle with a refraction 
index n > 1, and a massive mode corresponding to a massive particle with mass of the 
order of the exchanged scalar/pseudoscalar one. By taking into account the effect of the 
scalar/pseudoscalar width F, we will show that the massive propagation mode is a solution 
of the dispersion relations only when the size of F would be smaller or comparable to the 
effect of the mixing term, or in other words the massive solution is allowed only when the 
external magnetic or electric field is above some particular critical value. 

Finally, we would like to stress that the scalar/pseudoscalar coupling with two photons 
can induce a sizeable contribution to the photon splitting process 7 —^ 77 in external 
magnetic or electric fields, depending on the magnitude of the external field and the 
photon energy involved. However, although suppressed by the corresponding width of the 
scalar or pseudoscalar field, the corresponding photon absorption coefficient may turn out 
to be significant due to a strong resonant phenomenon when the external field is in the 
vicinity of the critical value. 

These results suggest a new class of experiments for searching indirect effects of the 



scalar/pseudoscalar coupling to two photons, like for instance by searching for the photon 
splitting 7 — i> 77 in constant magnetic fields. Notice that this process does not have any 
background since the analogous effect in QED is strongly suppressed for laboratory mag- 
netic fields and photon energies in the optical or X-ray frequency range [21]. Therefore, 
any observation of photon splitting signal in laboratory would be a clear indication of new 
physics effects. 

Returning to the scalar coupling with two photon, a well-known example in the stan- 
dard model is provided by the Higgs boson [2]. Unfortunately, in this case a too large 
critical magnetic (or electric) external field would be required for photon energies of the 
order of TeV, with no practical laboratory applications. Nevertheless, we stress that size- 
able effects on the photon absorption coefficient, induced by the Higgs boson, might be 
possible in astrophysical context, where sources of strong magnetic fields and/or very high 
energy photons can be found. For example in the core of a neutron star magnetic field 
can be as large as 10^ Tesla [5], although in that case plasma effects should be taken into 
account. 

The plan of the paper is the following. In section 2, we present the results for the 
effective photon propagator in the presence of a static and homogeneous external magnetic 
field. In section 3, we provide the solutions to the photon dispersion relations, while in 
section 4 we will show the results of the photon absorption coefficient. In section 5 the 
phenomenological implications of these results will be analyzed and numerical results will 
be provided. Finally, in section 6 we will present our conclusions. In appendix A we 
report the exact solutions of the dispersion relations, while in appendix B we provide the 
calculation of the imaginary part of scalar(pseudoscalar) self-energy. 

2 Effective photon propagator 

Let us consider first the case of a scalar field and an external constant and homogeneous 
magnetic field B assumed to be extended in all space. The Lagrangian containing the 
pure radiation field and other dynamical fields, can be obtained by decomposing the 
electromagnetic field strength F^^ in two parts 

F,, = F,M) + FiS^''' (3) 

where F^y{A) = d^A^ — di^A^ contains the usual photon radiation field A^, while the 
F^^^^ includes the corresponding term induced by the external magnetic field. After the 
shift in Eq.(3), the relevant Lagrangian would be given by 

L = -i F,M)F^'^{A) + ^d.^sd^^s - ^rnVs - ^^F,^{A)F^''{A)^s + ^cxt(B) (4) 

where m is the mass of the scalar field and Lcxt(B) is given by 

i^ext(B) = _iF^-t)^(ext)^. _ l^^ext)^^.(^) _ _i_^M^(^)^M)^^ 



Notice that the last term in Eq.(5), in the case in which the external field is constant, is a 
tadpole term for the scalar field ips- It is easy to see that this tadpole has no physical effect 
since it can be eliminated by making the following constant shift ips ^ fs + ^ where 6 = 
_^M)i?{ext)Mi^/(^4^^^2^)_ r|nj^g g^^^^ ^g^^ __|_p^^(^A)F^"'{A) can be then re-absorbed 

in the photon wave function (A^) renormalization A^ -^ Z^^'^ A^ with Z = (l + — j 
There are no other effects of the shifting, apart from adding an extra constant term to 
the Lagrangian.^ 

Finally, the interaction Lagrangian which is relevant to our problem is contained in 
the third term in Eq.(5), in particular 

4^„t(B) = -^ B ■ (v. A l(x)) cpsix) (6) 

where the symbol (a A b) indicates the standard vectorial product, V^: = d/dxi, i = 
(1,2,3) and (ps{x) and A^ are the scalar and photon fields, respectively, with A{x) = 
Ai{x), i = (1,2,3). In momentum space, the corresponding Feynman rule associated 
with the interaction in Eq.(6) is given by 

OO^^I (k) 



where Eimn is the antisymmetric tensor in 3 dimensions, k^ is the (m = 1,2,3) spatial 
component of the photon A;^ = {uj, k) entering in the vertex with polarization vector ei{k). 
Notice that at the vertex the energy-momentum is conserved, and that the external field 
B does not carry any four-momentum (p^ = 0), since it is assumed to be space-time 
independent. 

We are interested in analyzing the modification of the photon dispersion relations 
induced by the interactions in Eqs.(l),(2), in the presence of a constant and homogeneous 
external field. To a first inspection of the leading order effects to the photon self-energy, 
these contributions split in two separate class of diagrams: 

• the one-loop diagrams with no external field insertions induced by the operator 
F^,{A)F^^''{A)ips; 

• the tree-level diagrams with two external field insertions induced by the mixing 
operator F'^^F^^''{A)lps , see Fig.l. 



^Regarding the mixed term F^J^ 'F'^'^{A), we will not consider it in the analysis since it will affect 
only the exchange of photons with zero energy and momentum, and so would not have any effect on the 
dispersion relation of photons with frequency w > 0. 






Figure 1: The effective ptioton propagator obtained by summing up the full series of diagrams 
with external field insertions (curly lines with a cross) carrying zero momentum. The big bubble 
in the middle of the dashed lines indicate the exact summation of the scalar/pseudoscalar self- 
energy. 

The contribution of the first class of diagrams described above, since there are no external 
field insertions, cannot modify the dispersion relations. This last effect can only indeed 
be achieved by summing up the contribution of the second class of self-energy diagrams. 

For this reason, in the effective photon propagator, obtained by summing up the self- 
energy diagrams to all orders, we do not include the contribution coming from the first 
class of corrections induced by the operator Ffj^^{A)F^'^ {A)ips to the photon self-energy. 
Their effect, as well as the analogous one of standard QED vacuum polarization diagrams, 
can be indeed re-absorbed in the photon-wave-function renormalization. However, there 
are always higher order contributions coming from the mixing corrections induced by the 
F^^{A)F'^'^ {A)(fs operator and the vertex interaction in Eq(6). These higher loop effects 
should be considered as next to leading order corrections to the self-energy contribution 
induced by the term in Eq(6) and we will neglect them in our analysis. 

Now we consider the contributions to the effective photon propagator obtained by 
summing up the full series of diagrams as in Fig.l. In our calculation we include also 
the renormalized self-energy effects in the scalar propagator, summed up to all orders, as 
indicated in Fig.l with a bubble inside the scalar propagator (dashed line). Notice that 
the leading contributions to the scalar self-energy are provided by the one-loop diagrams 
in which two photons lines are circulating in the loop. As shown later on, these class of 
diagrams will be crucial for our analysis since they induce a (finite) imaginary part in the 
scalar self-energy as well as in the effective photon propagator. 

Let us choose a gauge where the free photon propagator Dij{k) in momentum space 
has only spatial components 

DM = p^^.. , (7) 

and zero for all the other ones. From now on, with the symbol k"^ we indicate k^k'^, 
where k^ = u and ki = {k)i are the corresponding energy and 3-momentum, respectively, 
associated to the photon propagator, and 6ij is the unit matrix in 3-dimensional space. 
In this gauge, the result for the effective photon propagator P-j ' (A;,B), induced by the 



mixing with a scalar field in the presence of an external magnetic field B, is the following 

where the tensor functions T/- (B) and -R^-- (B) are symmetric tensors in the i, j indices, 
and T/j -^(B) is a projector for k and B vectors, that is satisfying the conditions kiT^^ (B) = 
BiT//\B) = 0. More explicitly, 

(1)^ \B\'^kikj -k-B {kiBj + kjBi) + \k\'^BiBj 

^i (B) = 5ij ^-^ ' ^ ^.2 

IfcPlBP- [k-B] 



^(1) ,^, |BpA;,fc,- -k-B {hB, + k,B,) + |A;|2B,B, 

^ij (Bj = ^ z .^ ^.2 ' \^) 

IfcPlBP- (A;-B^ 



and the expression for scalar part of photon self-energy 11*^^^ (A;, B) induced by the external 
magnetic field B is given by 

nw, -X 1 IA^PIBP sin^^ , , 

where ms and T^sik"^) are the renormalized mass and self-energy of the exchanged scalar 
boson, and 9 is the angle between the direction of the photon momentum k and the 
external magnetic field B. Notice that the first term on the right hand side of Eq.(8), 
which includes the contribution of longitudinal polarizations, depends on the gauge choice 
in Eq.(7), whereas the second term in Eq.(8) is gauge independent. The latter property 
comes from the transversality of T^ • (B) tensor, namely kiT/^, (B) = 0. 

The physical interpretation of the T/^- (B) and Rl^'iB) is clear. The Rl^ {B) represents 
the contribution of both longitudinal and transverse photon polarization components, 
while T^j (B) represents the contribution of the transverse polarization orthogonal to 
both B and k. This picture is largely simplified by choosing the momentum k direction 
orthogonal to the external field that is A; ■ B = 0. Then 



' ' |A;|2 |B|2 '' ^ ' \k\^ |B|2 



As expected from the structure of the interaction vertex, the external magnetic field 
modifies the photon dispersion relations only for the transverse polarization component 
which is orthogonal to the plane generated by external magnetic field direction and the 
photon momentum. This modification is contained in the denominator of the second term 
of Eq.(8) due to the self-energy correction Ii'^^\k^ B) in Eq.(lO). 

Now we consider the case of a scalar field coupled to an external electric field E. By 
retaining only the linear terms in the external field, the relevant interaction Lagrangian 



is given by 



LLi^) = -^^^-(§-/i^))vsix). (12) 



The corresponding Feynman rule is in this case given by 
(k) I (k) 

Ai ;t^^.^ (p^ A^ ' 



where Ej is the i-th component of the external electric field E and u is the corresponding 
energy of the photon whose four-momentum is entering the vertex. Notice that the 
expression of the Lagrangian in Eq.(12) depends on the adopted choice of temporal gauge 
Aq = 0. However, it is worth stressing that the photon dispersion relations, being an 
observable physical effect, will not depend on the gauge choice. 

In analogy with the external magnetic field, the effective photon propagator induced by 
an external electric field can be straightforwardly derived by summing up the self-energy 
diagrams to all orders as above. By using the same gauge condition as in Eq.(7) for the 
free photon propagator, the result for the effective photon propagator is the following 

'' ^ ' k^ + ie '' ^ ' A:2-n(2)(A;,E) ' ^ ' ^ ' 

where T/^ (E) and /^^^(E) are given by 

T,f(E) = S|, flg)(E)=%-S|. (14) 

The tensor functions Tj^ (E) and R\j {E) represent the contributions of the photon po- 
larization parallel and orthogonal to the external electric field E, respectively. The scalar 
part of the photon self-energy in the second term of Eq.(13) is given by 

In the case of an angle 9 ^ tt/2 formed by the momentum k of the incident wave and the 
direction of the external electric field E, the formula in Eq.(15) will be modified as 

n(2) ik, E = ^--r — ^^ —r- . (16) 

and the corresponding photon polarization will be the one parallel to the plane formed 
by the momentum k and the external field E. 



Contrary to the case of an external magnetic field, the external electric field modi- 
fies the photon dispersion relations only for the photon polarization component which is 
parallel to the plane formed by E and k, while the orthogonal component remains unaf- 
fected. The modification of dispersion relations is contained in the denominator of the 
second term of Eq.(13) due to the self-energy correction n*^^)(A;, E) in Eq.(16). 

Analogous results for a pseudoscalar field minimally coupled to two photons as in 
Eq.(2) can be obtained in a straightforward way from the above results. In particular, for 
a pseudoscalar field the relevant interaction Lagrangian, linear in the external field, these 
are given by 

4nt(E) = -^ E . (v. A Aix)) Mx) , (17) 
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^h.t(B) = -— B. -A(x)Up(x). (18) 



By comparing the Lagrangians in Eqs.(6), (12) with the corresponding ones in Eqs.(17), 
(18) it is easy to derive the effective photon propagator by using the previous results for the 
scalar case. In particular, for the effective photon propagator induced by a pseudoscalar 
field in the presence of external magnetic or electric field, one obtains 

P^T'^\k,B) = yy^— i?!f(B) + ^-^T,f(B), (19) 

^ k^ + te '^^ ^ F-n(2)(A;,B) *^ ^ ^' ^ ^ 

P,f'^)(A;,E) = ^_i?«(E) + '- ^T^I\e), (20) 

' k^ + ie '' ^ ' A;2-n(i)(A;,E) ' 

provided that inside the functions t//^(E), i?Jj^(E), t//^(B), i?Jf (B), n(2)(A;,B), and 
n*^^)(A;,E), appearing in Eqs.(19),(20), the following substitutions A^ -^ Ap, ms -^ mp, 
and 'Es{k'^) -^ 'Ep{k'^) for the corresponding pseudoscalar quantities are implemented. 

3 The photon dispersion 

In order to derive the new dispersion relations for the propagating photon induced by the 
external field, we first look at the zeros connected to the inverse of the real part of photon 
propagators in the last terms of Eqs.(8), (13) and (19), (20). In particular, this consists 
in solving the following mass-gap equations^ 

^Even if our mass- gap equations allow for a massive solution of the type k^ = M^, this should not 
be connected to the presence of a phase transition like in the BCS theory of superconductivity or in 
the Nambu Jona-Lasinio model. A mass term is already present in our Lagrangian and it is given by 
the scalar/pseudoscalar mass term. The effective photon propagator can have a massive pole due to the 
presence of photon-scalar/pseudoscalar mixing term. Thus, the mass-gap equations above do not have 
the same meaning as in the dynamical mass generation mechanisms mentioned above. 



Scalar 



External B -^ P - Re [nW(A;, B)] = 0, 

External E -^ k"^ -Re[U^^\k,E)] = 0, (21) 



Pseudoscalar 



External B -^ P -Re[U'^^\k,B)] = 0, 

External E -^ A;^ - Re [n(^)(A;, E)] = 0. (22) 

The fact that there is the presence of an external field, does not guarantee anymore that 
the equations above have a unique solution at /c^ = 0, as expected instead in the case of 
vanishing external sources. 

However, we will see that not all the solutions of Eq.(22) will correspond to physical 
ones, and there will appear a spurious solution. In particular, in order to link the solutions 
to the physical spectrum, the following conditions must be fulfilled: 

• the square of the photon energy u must be real and positive, namely u"^ > 0; 

• the residue at the pole of the propagator must be positive, in order to ensure that 
the corresponding quantum state has a positive norm. 

We will show that the above equations admit at most two acceptable solutions which 
satisfy the above criteria. This can also be seen by switching off the imaginary part of 
the scalar /pseudoscalar self-energy. In this case the equations above become quadratic in 
/c^ and only two solutions are generated, corresponding to the classical ones [16]. One of 
these solutions corresponds to a light-like mode {u = \k\/n) with a refraction index n > 1, 
while the other one can be associated to a massive mode (cj^ ^ \k\'^ + M^) with a mass 
M very close to the mass of the exchanged scalar or pseudoscalar particle. Indeed, the 
total number of degrees of freedom must be conserved. In the original Lagrangian there 
are three degrees of freedom corresponding to particles which are on-shell, the two photon 
components and the scalar/pseudoscalar field. The number of physical poles are indeed 
three, the two mentioned above and the pole at /c^ = related to the photon polarization 
which does not interact with the external field. 

As we will show later on, the inclusion of the quantum corrections, which are absorbed 
in the scalar/pseudoscalar self-energy, cannot be neglected due to the presence of an 
imaginary part in the scalar/pseudoscalar self-energy S(A;^) induced by the interaction 
in Eqs.(l) and (2). Then, the massive dispersion mode would be allowed provided that 
the external field is above some particular critical value, which depends on the mass and 
decay width of the exchanged scalar/pseudoscalar particle and photon energy. 

Now we will provide the solutions to the mass-gap equations in Eqs.(21), (22). The 
exact results are given in the appendix A, while below we will report the approximated 
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expressions obtained by expanding the exact solutions in powers of a small parameter 
contained in the photon self-energy 11*^^'^) (A, A;). In order to simplify the notation, it is 
convenient to compactify the set of Eqs.(21) and (22) as follows 

A;=^ - Re [n(A;, A)] = (23) 

where n(A;, A) generically stands for the effective photon self-energy and it is given by 

nftA) = ^,_,/_^(^,) ^ (24) 

Above, m and S(A;^) indicate the renormalized mass and self-energy of the exchanged 
scalar or pseudoscalar particle, respectively. In Eq.(24), the term A absorbs the contri- 
bution of the external field interaction and it is given by A = A*^^)(B) and A = A*^^^(E) 
for a scalar field in the presence of an external magnetic and electric fields, where 



lA^PlBpsin^^ .2UA tu^lEpsin^^ 



a^^hb) = '"";:' r , a^^he)= " z7 (25) 



and similarly, A = A*^^)(E) and A = A*^^)(B) for the pseudoscalar case. 

By separating the real and imaginary part in Eq.(24), the Eq.(23) can be written as 

^^, A(P-m^-Re[S(P)]) _^ 

(P-m2-Re[E(P)])' + (Im [S(P)])' ' ^ ^ 

When /c^ = m^, the term m'^ — Re[S(m^)] corresponds to the renormalized scalar or 
pseudoscalar mass, including radiative corrections. However, it is convenient to adopt the 
so-called on-shell mass renormalization scheme where Re[S(A;^)] 1^2=^2 = 0. Then, from 
now on, m^ will stand for the renormalized mass. On the other hand, Im[S(A;^)] is finite 
at 1-loop, being connected to the tree level decay in two massless photons. In particular, 
see Appendix B for more details, for scalar (pseudoscalar) interaction we have 



647rA 



M^k^)] = -7,h^ ■ (27) 



S(P) 

We stress here that the threshold where the imaginary part of the scalar /pseudoscalar 
self-energy Im[E(A;^)] is different from zero starts at k"^ = 0, due to the fact that a 
scalar/psedusocalar particle with two photon interaction can always decay in two massless 
photons, regardless of the size of its mass. We will see that the effects of the external field 
corrections, will shift this threshold above k"^ > due to the modification of the photon 
refraction index, see section 4.2 for more details. 

In the physical scenarios that we are considering here, there is a characteristic (small) 
dimensionless parameter which is given by A/m'^. Its smallness is due to the fact that the 
A ~ 1/A|p where the effective scale As^p is much larger than any other characteristic 
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energy scale present in our problem, namely A|p ^ I |B|, |E|,cj^,m^ i. This observation 
will simplify the analysis regarding the solutions of Eq.(26). An approximate solution 
can indeed be easily found by using the Taylor expansion in terms of A/m'^. By looking 
at Eq.(26) one can see that there is always a solution around the massless pole /c^ ~ 0, 
which can be parametrized as k"^ = Eq, with Sq/tii? of order A/m^. This corresponds to 
the dispersion mode where the effects of external field modifies the refraction index of the 
vacuum, as for instance in the case of QED vacuum polarization in the presence of an 
external magnetic field [21]. However, as discussed before, Eq.(26) admits also a massive 
solution provided that A/(Im[S(m^)])^ ^1. As we will show later on, this inequality 
will imply the existence of a critical value for the external field. 

The Eq.(26), being a cubic equation in A;^, admits three solutions, namely /c^ = 
Mq, M|.. The general expressions are given in the appendix A. However, approximated 
solutions can be easily found by using the following antsatz 

k^ = So (28) 

fc2 = m'^ + Sm, (29) 

By substituting these expressions in Eq.(26) and neglecting higher order terms 0{eQ/m'^) 
and 0{e'^/m'^), the algebraic equations in k'^ can be split into a linear and quadratic ones 
for £0 and Em, respectively. We get the following results 

A 
A 



where 



Mo = 2 



Ml = ^' + ^(l±\/l-e), (30) 



^-^^' (31) 

and r is the scalar (or pseudoscalar) width given by Im[S(m^)] = —mT. As shown in 
the appendix A, these results can also be easily obtained from the exact expressions by 
retaining the leading order terms in the I\jvr& expansion. In the value of F we absorb 
all the effects of the external fields corrections on the scalar /pseudoscalar decay width. 
However, it is reasonable to assume that these corrections are sub-leading for our problem 
and for practical purposes we neglect them and retain only the tree-level contributions 
to the scalar /pseudoscalar width induced by the interactions in Eqs.(l),(2). In order to 
have real solutions for m^, the condition ^ > 1 must be satisfied. We will return on this 
point after discussing the dispersion relations. 

As previously discussed, the lightest mode (Mq) corresponds to the propagation of a 
photon in a medium with a refraction index n given by 

rsw-f^N . |BP sin^6' ,„^, 
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n 



(s) 



(E) 



1 + 



|E|2 sin^^ 
2m2|A5|2 



(33) 



where n is usually defined as n = \k\/iu. Above, n^ '(B) (w '(E)) stands for the refraction 
index in magnetic (electric) external fields with scalar interactions. Analogous results are 
obtained for the pseudoscalar case, with n*-^-'(B) = n*-'^-'(E) and n'^^^'E) = n^^\B) , where 
the substitution of A5 — > Ap in n^^\B) and n^^^E) is understood. The next leading order 
corrections in A/m^ to the dispersion relations gives the non linear dependence of the 
refraction index with the photon energy. From now on in our notation, whenever the 
energy {u) dependence in the refraction index n{uj) is not explicitly shown, it means that 
it corresponds to n(c<j = 0). 

Regarding the other massive solutions k"^ = M^ of Eq.(26), only one is physically 
acceptable and it will correspond to the pole k"^ = M^. This can be easily seen by taking 
the limit F ^ or analogously ^ — > 0, where one should recover the classical (F = 0) 
solution /c^ = m^ + ^ [16]. This massive pole should correspond to the propagation 
of a particle of mass m = M_|_ provided that the square of the mass term M^ is a real 
and positive quantity. At this point it is convenient to define the following dimensionless 
parameters xe,b by xb = ^mZA^ ^'^'^ ^e = 2m?A^' ^'^'^ problem xe,b ^ 1 and 
one can simplify the dispersion relations above by expanding them around xe,b = 0. In 
particular, for the scalar case, by using the approximate solutions in Eq.(30), we have for 
the massive mode in the presence of external electric and magnetic fields 



External B 



up ^ 



k\^{l+XB\l 



4F2m6A4 



\ |B|4sin^^|A?|4^ 



+ m 



(34) 



External E 



uj"^ ^ 



m 



1 + Xe 



V 



\ 



V^m? 



X, 



{\k 



2 _|_ ^ 



-\] 



J ) 



(35) 



where terms of order 0{x% ^) have been neglected. Analogous results for the pseudoscalar 
case are obtained by using for the external magnetic and electric fields the expressions in 
Eq.(35) and Eq.(34), respectively, and by substituting A^ -^ Ap. 

One can easily see that, in order to have a real solution for o;^, see Eq.(30), the 
following condition must be satisfied 



A > 2m^F. 



(36) 



Notice that in the classical approximation [16], no condition is required for all the pa- 
rameters in the Lagrangian in order to have a massive pole in the spectrum. Indeed, 
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if one sets F -^ 0, the condition (36) is always satisfied. However, we emphasize that 
the approximation of neglecting the width F in the axion physics, as usually done in all 
previous studies, is well justified due to the very small axion mass. Indeed, for character- 
istic magnetic fields of the order of Tesla, axion masses of the order of meV, and photons 
wave lengths in the optical region, the condition (36) is always satisfied, provided that 
the contribution to the total width is dominated by the axion decay in two photons. 

The physical meaning of Eq.(36) can be roughly understood as follows. Let us consider 
the case in which the width F is very large in comparison with the mass scale set by the 
quantity A, in such a way that the relation (36) cannot be satisfied. In this case, the 
fact that the scalar/pseudoscalar particle can decay too fast does not allow the massive 
photon mode to be formed and hence cannot coherently propagate. 

Given a particular value for the photon energy u, the inequality above leaves to the 
following critical values for the magnetic and electric external fields for a scalar interaction 

mi\^\As\ I 777 

|B| > Bent^^ \-^Jl + R{5) 

\ksm.9\ 

mi V2r|A. 



a; I sin 6' 
and analogously for pseudoscalar one 



E| > E,,it = "^ :" ::''V i + i?((5) (37) 



m2V2F|Ap| r 



|E| > Eent = ^ '-^Jl + R{5) 

fcsin^ 



|B| > Be.it ^ ^';^^'f^' yrTi?(^. (38) 



where \k\ ~ yuP^^^rri? + 0{xb,e)- The expression for the function R{x) ~ 0{x), con- 
taining the higher-order corrections in powers oi S = T/m, is reported in the appendix A. 
Analogously, for a fixed value of the external field, the conditions above can be read as 
the minimum photon energy ui necessary in order to generate a massive mode. We have 
explicitly checked that the approximated solutions in Eq.(30) are obtained from the exact 
ones by retaining only the leading contribution in the weak external field expansion. See 
Appendix A for more details. 

A remarkable aspect of the results in Eqs.(37), (38) is that in the case in which the 
scalar or pseudoscalar field has only the decay mode in two photons, the critical value 
for the external field does not depend on the coupling A of the effective interaction. This 
can be easily checked by noticing that the corresponding width F in that case would be 
proportional to F oc rrt" / A?, while A oc 1/A^. 

In order to identify the massive poles in the effective propagator as physical quantum 
states, one has to check that the corresponding residue at the pole of the propagator is 
positive. Indeed, the residue at the pole {Z) is connected by unitarity to the norm of the 
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quantum state excited from the vacuum. In particular, for a generic solution k^ = M^ of 
the mass gap equation, one gets 

Z= lim (e-M^)(P-ReU{k,A)y\ (39) 

or equivalently 



d_ 




^=(l-7^Ren(A;,A)_J . (40) 



where Z^/^ is the wave function renormalization constant. 

By using the expression in Eq.(40), we find the following result for the residue at the 
poles k"^ = Mq and k"^ = M]_, respectively Zq and Z± 

Zo = 1-4' (41) 

Z± = ±^- ,- ^ , (42) 

where we have neglected terms of the order 0{^). As we can see from these results, two 
physical solutions are allowed, corresponding to the positive norm states oi Zq > and 
Z^ > 0, connected respectively to the poles k"^ = Mq and k'^ = M^, while k'^ = M^ is an 
unphysical one, being associated to a ghost (Z_ < 0). 

Notice that when the external field is far above the critical value, .^ <^ 1, but still in 
the weak coupling regime i.e. A/m^ <^ 1, then the term Z+ ~ (9(A/m^) <^ 1. This 
means that in this case the probability to induce the massive mode from the vacuum, is 
very suppressed. However, when we approach the critical value C ^ 1 the constant Z^ 
increases and tends to infinity at .^ = 1. 

This behavior can be easily understood if we look at the inverse of the propagator 
n~^(A;^) near the pole k"^ ~ M^, in particular 

U-\k^)\k2^M^ = (P - M^)Z-^ + (P - M2)2ReS"(M2)|fc2=M2 + • • • + ilmi:{M^) (43) 

where the dots stand for higher order terms in {k^ — M'^)/m'^ expansion, and S (M^) 
indicates the second derivative of S(A;^) with respect to A;^ and evaluated at A;^ = M^. 
Notice that .^ = 1 corresponds to two degenerate massive solutions, and the propagator 
generates a double pole in /c^ = M^ = M^. This is not a matter of the approximation 
adopted, since the same conclusions are obtained in the exact case (see Appendix A for 
more details). Clearly, since Z~^ is related to the coefficient of the single pole, Z~^ 
vanishes at .^ ^ 1 or analogously Z ^ oo. Moreover, since the Z factor is related by 
unitarity to the probability of creating the corresponding quantum state from the vacuum, 
when Z > 1 the unitarity is spoiled. Therefore, the requirement of unitarity sets a region 
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of validity of our calculations when the external field is close to its critical value. By 
imposing Z^ < 1 we obtain 

< e < 1 - — + 0(— ) . (44) 

In conclusion, the probability that the massive mode is excited from the vacuum increases 
as the external field approaches (from above) the critical value, but clearly vanishes just 
below the critical value. 



4 The photon absorption coefficient 

Due to the presence of an external constant electric or magnetic field, the interactions 
in Eqs.(l), (2) can induce a non- vanishing probability for the photon splitting amplitude 
7 -^ '-ff. This phenomenon is well known in the framework of QED [21], as well as the 
analogous one of electron-pair creation 7 — ^ e+e^ in constant and homogeneous magnetic 
field [22,23]. Nevertheless, there are no studies so far concerning the analogous effect of 
photon conversion process 7 ^ 77 induced by the scalar and pseudoscalar interactions.^ 
Measurement of a photon splitting in constant magnetic field is a challenge for experiment, 
although high-energy photon splitting in atomic fields has been recently observed [25]. 

For this kind of problem it is appropriate to express the conversion probability in 
terms of a photon attenuation coefficient, usually called absorption coefficient (dy) [23]. 
In particular, the number A^^^ of photon splitting events created by a photon crossing an 
EM background field for a path length L are given by 

N^^ = N^{1- exp [-dy L]) ~ N^d^ L (45) 

where A^^ is the total number of photons entering the background EM field. Notice that 
the last relation is a good approximation only when d.yL <^ 1. 

Due to the photon interactions in Eqs.(l), (2), a new kind of contribution to the 
photon absorption coefficient d^ is expected with respect to the corresponding vacuum 
polarization effects in QED [21-23]. Now, before entering into the analysis of the new 
physics contributions, let us start by recalling the known results on the photon propagating 
in a constant magnetic field [21]. 

In QED, the matrix element for the 7 — ;> 77 process can be calculated in perturbation 
theory by using the Euler-Heisenberg effective Lagrangian [26] arising from integrating 
out electrons at 1-loop. After expanding the electromagnetic field around the constant 
background field, new dispersions relations are obtained for the photon propagating in the 

•^However, there is an analogous study analyzing the effects of external field on photon and axion (a) 
decays 7 — ^ 07 and a -^ 77 in the framework of a very light axion [24]. We stress that the analysis 
and results contained in [24] are quite different from the ones presented in our work and there is not any 
significant overlap. 
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external field. A refraction index (n(A) > 1) associated to the photon polarization £^(A), 
where A parametrizes the two polarizations parallel and perpendicular to the magnetic 
field direction, arises. Even by neglecting the dispersion effects, the photon splitting 
process 7 — > 77 is kinematically allowed, provided that all particle momenta in the 
reaction are proportional to a unique momentum k^, satisfying k"^ = 0. However, as 
shown in Ref.[21], at least three external field insertions would be necessary in order to 
have a non- vanishing matrix element. This is a simple consequence of the Lorentz and 
gauge invariant structure of the Euler-Heisenberg Lagrangian and due to the fact that 
there is only one light-like four-momentum in the process. The magnitude of the resulting 
photon splitting absorption coefficient would then given by [21]: 

^r-o.i(^)°(i^)\,„-. (4a) 

where B^^^ = ^^, with e the unity of electric charge, is the critical magnetic field for 
the photon-pair creation [22, 23] and rrie is the electron mass. This formula holds only 
for 6 = 7r/2 and in the weak- field limit eB/ml <^ 1. Then, one can see that even for 
strong laboratory magnetic fields (of the order of Tesla) the absorption coefficient d^ 
would be very small, due to the fact that B^^^ ~ 10^° Tesla. Only in astrophysical 
context this effect can become relevant. In particular, with pulsar magnetic fields of the 
order of B^^^ and uj ^ me, there are many photon splitting absorption lengths in a 
characteristic pulsar distance of lO^cm. By taking into account the dispersive effects, the 
energy-momentum conservation is modified. The momenta of the two final photons will 
not be parallel anymore and differ from each other for small angles. Due to the dispersive 
effects, the matrix element of the photon splitting can be induced by one external field 
insertion. However, in this case the corresponding d^ would turn out very suppressed by 
kinematical factors being proportional to the small opening angle of the final photons. 
[21]. 

Returning to our case, a practical way to calculate the absorption coefficient from the 
photon self-energy in Eq. (24) is by making use of the optical theorem which connects d^ 
to the imaginary part of the self-energy evaluated at the corresponding poles of the prop- 
agator. As can be seen from Eq.(lO), the imaginary part of photon self-energy ImH(/c, A) 
is then proportional to the imaginary part of the scalar/pseudoscalar self-energy, namely 
Imn(A;,A) ocImS(A;2). 

Now, if we do not take into account the effects of vacuum polarization in background 
EM field, and consider the photon purely massless, the kinematics of the reaction 7 — » 77 
would not forbid this process provided that all the momenta are proportional to a unique 
light-like momentum k'^ = 0. However, as in the analogous case of QED, due to the fact 
that there is only one independent light-like four momentum and due to the antisymmetric 
property of the electromagnetic field strength F^,y, the minimum number of external field 
insertions is equal to three in order to have a non vanishing effect in the matrix element. 
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This is a general result, as proved by Adler [21], and does not depend on the structure 
of the interaction, but on the gauge and Lorentz invariant property of the effective La- 
grangian for the photon obtained after integrating out the other degrees of freedom. In 
the case of weak magnetic fields, this would lead into a very strong suppression. As we 
will show later on, if dispersion effects are taken into account, (/c^ = Mq), the process 
could proceed by means of one external field insertion only, but suppressed by the small 
angle induced by the dispersions effects. In the next two sub-sections we will analyze 
the photon splitting phenomenon induced by the massive and massless dispersion modes 
respectively. 

4.1 Massive dispersion mode 

As shown in the previous section, the massive mode of the photon, with mass of the order 
of the scalar/pseudoscalar one {m), is allowed by dispersion relations when the critical 
conditions in Eqs.(37), (38) are satisfied. Then, the m,assive mode could easily decay 
in two lighter photons, or in any other kinematically allowed final states / coupled to 
the intermediate scalar/pseudoscalar particles. In this case, the corresponding width will 
be proportional to the imaginary part of self-energy evaluated on the massive pole of 
the photon propagator (A;^ = M_^). In conclusion, the photon could get a non- vanishing 
decay width (F^), provided that the critical conditions in Eqs.(37), (38) are satisfied and 
u > m^, where m^ = J M\. 

One can formally define a width (F^) associated to the m,assive mode by making use 
of the optical theorem. In particular, 

m-,F^ = -Z+Im[n(A;,A)] ^^_ ^9{ij - m^)9{B -B,,,t) , (47) 

where 6(x) is the standard ^-function defined as 6(x) = 1 for x > and 6(x) = for 



a; < 0, and m^ = Jm? + 2^(1 + VI ~ 0- The constant Z+ appearing in Eq.(47) is the 
corresponding residue at the pole, given in Eq.(40). As shown later, the same result for 
F^ in Eq.(47) can be re-obtained by starting from the matrix element of the transition 
7* -^ 77, where 7* represents the massive mode of the photon. For the imaginary part 
of photon self-energy evaluated at the massive pole k'^ = mi we have 



Im [n(A;, A)] |,2=^2 = -m^A 1 . (48) 



Then, the F^ is given by 



^-TT^i-^^^^^^l' 
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where the last term in parenthesis conies from the first order expansion in A/iri^ of the 
m^ when expressed in terms of the axion mass m. The absorption coefficient d^ entering 
in Eq.(45), is related to the F^ in Eq.(47) by 

d^ = 1}lT^. (50) 

Finally, by inserting the results of Eqs. (42), (47), and (48) into Eq. (50), we obtain 

where we neglected terms of the order OiV"^ /m?). We stress that to be more precise 
the scalar/pseudoscalar width F appearing in all the calculations above should be the one 
evaluated on the massive pole m^, that is F(r7i^), which differs from the width F evaluated 
on the scalar/pseudoscalar mass-shell mass m by small terms of order 0{T /m). The same 
argument applies for the mass m appearing in the expression mV which should be m-y 
rather than m. 

As explained at the end of the previous section, the restriction on the upper limit 
^max _ ]^ _ i_ comes from the requirement of unitarity Z^ < 1. Indeed, the validity of the 
results in Eq.(51) is based on the optical theorem, which holds only under the hypothesis 
of unitarity. Notice that, for the imaginary part of photon self-energy evaluated on the 
other pole /c^ = M^, the corresponding absorption coefficient would have been negative 
due to the fact that Z_ < 0, pointing out the presence of an unphysical solution. 

These results can be easily re-obtained by starting from the standard formula for the 
absorption coefficient [21] of photon splitting 7(1-0') -^ '-y{uji)'^{uj2) 

where \M('j -^ 77) P is the square modulus summed over final state polarization of the 
corresponding amplitude. The extra factor 1/2 in front, takes into account for the identical 
final states of two photons. The amplitude at the leading order can be easily obtained by 
evaluating the following Feynman diagram 

<Y(ki) 
V 
y(^2) (53) 

where the dark and light curly lines represent the massive and massless photon propa- 
gation modes of the photon, and the bubble stands for the photon-scalar/pseudoscalar 
vertex induced by the external field, and the dashed line represents the scalar/pseudoscalar 
propagator. By taking into account the renormalization of the wave function Z_|_ in the 
massive mode, as given by Eq.(42), one obtains 
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where V^^{ki,k2) stands for the vertex of scalar/pseudoscalar in two photons and the 
sum is performed over the final photon polarizations. In particular, for the scalar (pseu- 
doscalar) interactions, V-y^{ki,k2) = V^^^\ki,k2), with 



1 
2X 



V^'^{k,,k2) = ^r^[{k,-k2)g^^ -k^,k^,]eUki)e2,{k2) 



V;^(A:i,A;2) = -^ (k^^k^ e^) e,,^p , (55) 

where Si{ki) and £2(^3) ^-^e the polarization vectors of the two final photons, and ki^2 
the corresponding four momenta. Then, by evaluating k"^ on the massive pole k"^ = m'i, 
Eq.(52) can be re- written as 

d^7 = r^ \f^ r2 2 I dH27T)'5\k -k,- k2) E \V,,{k,, A;2)p I , (56) 

where d<^ = 2b I (2^)320; (2w)^2uj ■ Notice that the integral over the phase space in the last 
term in parenthesis gives just the scalar/pseudoscalar width ^^r(m^) ~ — F. Finally, by 
integrating Eq.(56) and using the identity 

Z^A 1 



(m^ — m'^y + F^m^ ^1 — ,^ ' 



(57) 



one can easily recover Eq.(51). It is worth noticing that the denominator of the right-hand- 
side (r.h.s.) of Eq.(56), which is connected to the scalar (pseudoscalar) propagator, is going 
in resonance since m^ is very close to m. This effect partially removes the suppression 
given by the Z^ in the numerator, that is the probability to induce the massive mode. 

At this point it is fair to say that the results obtained for the absorption coefficient are 
based on the assumption that the EM background field is constant and homogeneous in all 
space. This is an approximation since for any practical experiment the external magnetic 
or electric field has finite extension. However, it is reasonable to believe that the effects 
of the boundary conditions can be neglected (as in our case) when the extension length L 
of the external EM background field where the photon is traveling satisfies the condition 
L ^ A^, where A^ ~ 1/m^ is the de Broglie wave length associated to the massive mode 
k'^ = m^. 

Now we analyze two particular limiting cases, the region of external fields very close 
to the critical value ,^ ~ 1, and the large external fields ^ ^ 1. It is easy to see that 
the maximum value of the width, compatible with unitarity, is obtained, as expected, 
near the resonant region ^ — * 1 where the Z^ -^ 00. By restricting ourself to the region 
of maximum value of Z^ allowed by unitarity, that is Z^ ~ 1 and corresponding to 
1^ ~ 1 — F^/m^, the maximum value of the absorption coefficient is 

m^ F 

' LU m 
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Figure 2: The ratio d^/dJ^^^, versus ^, for three values of (5 = F/m, with (5 = 10 ,10 ,10 
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In the opposite limit, of very large external fields ^ 
its minimum value given by 
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0, the absorption coefficient gets 



(59) 



which is independent of the external field. 

In Fig. 2 we plot the function d^ normalized to its maximum value d^^^ ^ versus ^, for 
some representative values of ratios Y jm. As can be seen from these results, the variation 
of the photon absorption coefficient as a function of ^ (or analogously as a function of the 
external field if all the other parameters are taken constants) is very steep only near the 
region close to the maximum value ^™^^ = 1 — F^/m^, when F/m ^ 1. Therefore, for 
very small values of Y jvn the average of d^ over the external field could provide a good 
estimation of the size of the effect, in particular 



Ci-Y 



c?^ d^ 



lYm 



UJ 



Y\ F2 



m 



m^ 



(60) 



However, notice that the average over the full available range of ^, (0 < ^ < 1), is ffnite 
and given by d^ = ^^-^(1 + 0{^)). The difference between this result and the one 
obtained in Eq.(60), with the restricted integration region is a quantity of order 0{Y/m). 
This difference should roughly give an idea of the theoretical uncertainty on the absorption 
coefficient when the average over the external field is considered. 



4.2 Massless dispersion mode 



Here we consider the possibility that the photon splitting in two photons is induced by 
the massless dispersion mode k'^ = Mq, see Eqs.(30) and (96). The transition amplitude 
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is given by the Feynman diagram as in Eq.(53), where now the initial photon dispersion 
corresponds to the massless mode. The evaluation of the absorption coefficient can be 
easily performed by using the standard formula in Eq.(52), as in the case of QED [21], 
provided that dispersion relations for the final photons are taken into account. At this 
point we should stress that if dispersion effects are neglected in the initial and/or final 
photons propagation modes (that is n = 1), then the corresponding absorption coeffi- 
cient vanishes due to the Lorentz covariant structure of the vertex V^^{ki, k2) in Eq.(55). 
Indeed, although the process of photon splitting with no-dispersion is not forbidden by 
kinematics provided that all momenta in the reaction are proportional to a unique mo- 
mentum satisfying /c^ = 0, the square amplitude X^poi 1^77(^1; ^2)!^ oc {ki • k2Y vanishes 
when k'l = kl = k"^ = 0. 

The standard formula for the absorption coefficient remains the same as in Eq.(52) if 
the dispersion effects are implemented in the delta function as [21] 

5'^{k — ki — k2) -^ 6{uj — Ui — uj2)5^{n{uj) cuk — ni{uji) cuiki — 722(^^2) ujk2) , (61) 

where u, uJi^2 and k, /ci 2 are the energy and modulus of the momenta of the corresponding 
waves respectively, k = k/\k\ is the unity vector associated to the momentum k, and 
generically n{uj) is the refraction index where the explicit dependence on the wave energy 
uj is shown. According to Ref. [21], the condition set by kinematics in order the reaction 
to proceed is 

S = ujini{uji) + LJ2n2{id2) - (uji + uj2)n{uji + LJ2) >0. (62) 

This condition can be easily understood by looking at the opening angle 612 between the 
two 3-momenta ki and k2 of the final photons, where [21] 

In order to have real values for 612, the condition S > must be satisfied. 

Let us start by analyzing the case of a photon crossing a constant and homogeneous 
magnetic field in the presence of a scalar interaction. In order to simplify the analysis 
we consider the case of a linearly polarized wave with momentum k perpendicular to the 
magnetic field direction. As a first approximation, we neglect the effects of the external 
field on the QED vacuum polarization. As we will show later on, this approximation 
can be well justified for typical mass and coupling A in the range of m ~ lO^^eV and 
A ~ lO^GeV, as suggested for instance by the results of PVLAS experiment [17], or in 
general for ymA < 5 x rrie ~ 2.5MeV, where rrie is the electron mass. 

From now on in this section we will use the following notation for the polarizations 
of the photons: namely 7|| (7±), corresponding to the polarization vector e\\(^±){k) parallel 
(perpendicular) to the plane formed by the external field direction B and the direction 



22 



of wave propagation k = k/\k\, where k ■ e{k) = 0. Moreover, n\\{u) and n±{uj) stand 
for the corresponding refraction indices. One can first notice that, in the case of scalar 
interaction, only the 7_l will acquire a refraction index n > 1, if the direction of k is 
perpendicular to B. 

Let us consider the dispersion relations for the polarization state 7^ with refraction 
index n±{uj) in the case of weak magnetic fields, or small frequencies, where A/m'^ <^ 1. 
At this point it is convenient to introduce an energy scale defined as /i = y/mA. As will be 
clear later on, this scale plays the role of an effective energy scale, as the analogous of the 
electron mass rrie in the QED photon splitting phenomenon. From Eq.(25), the term A 
(for 9 = 7i/2) is given by A = \k\^ \w~) where Bcr = /i^. ^ The solution of k"^ = Mq can 
be expanded in powers of photon energies u, provided that oj <^ jj,. By using the results 
reported in Appendix A, the dispersion relations at the next to leading order expansion 
in A/m^ give in this case 

-=«'(-(l;)"*S(£)><"(i)' 

The index n_[_{uj) is then easy to calculate and it is given by 

„,M^f:-f^v.^fi^vr" (65) 



where we neglected higher order terms in the -^ expansion. Three possible kinds of 

i>cr 

transitions for polarized states are allowed 

7± -^ 7± + 7± (66) 

7± -^ 7|| + 7|| (67) 

7± -^ 7± + 7|| • (68) 

In order to see if the reaction can proceed, one has to first check if the kinematical con- 
dition in Eq.(62) is satisfied. For this purpose it is useful to adopt the following approx- 
imation. Since the angles between the final and initial photon directions of propagation 
are very small, all the final momenta are almost aligned along the initial momentum. 
This implies that the modification of the refraction index n\\ corresponding to the parallel 
polarization wave 7|| of the final photon is much smaller than the corresponding one in 
the other polarization 7^, being suppressed by the small opening angle. In particular this 
means that ^||(0) — n^(0) < 0. Therefore, for the refraction indices of final photons we 
can approximate n\\ = 1, while keeping n±{uj) as the one in Eq.(65). By evaluating the 

''Note that the critical value of the external magnetic field Bcr associated to the massless mode is 
different from the critical magnetic field connected to the generation of massive mode, namely Bcrit, in 

Eqs.(37),(38). 
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expression S using the Taylor expansion up to the second order in the photon frequency 
[21], we get 

o 

S[7± ^ (7±) 1(7^)2] = -2^1(0) (^1^2 + ujiujfj > 

S[7± ^ (711)1(71)2] = (^1 + ^2) (n||(0) - nx(0)) < . 

S[7± -- (7||)i(7±)2] = cui(n||(0)-nx(0))<0, (69) 

where 

nl(0).|ln,MU„.-ij(|-)%0. (70) 

By requiring that S > 0, we see that the kinematical selection rule allows only for the 
polarized reaction 7_|_ — > 7±7±- We would like to stress that in QED the reaction 7_|_ — > 
7_l7_l is not allowed since n_|_(0) > 0, while in our case this is possible due to the fact that 
nl(0) < 0. 

Now we will provide an estimation of the absorption coefficient for the process 7_l -^ 
7_i_7_i_, by retaining the leading contributions in A/m'^ <^ 1 expansion. The square modulus 
of the amplitude is obtained by evaluating the diagram in Eq.(53) with scalar vertex 
interaction as in Eq.(55). By substituting Z^ -^ Zq and k"^ = Mq in Eq.(54) one gets 

|M(7x - 7±7x)P = [^y (i) ^ [(A;i ■ hM . e,) - (h ■ 6i)(4 ■ e^)f , (71) 

where higher order terms in A/m^ are neglected. In particular, we substituted the de- 
nominator in the right-hand-side of the Eq.(56) with the leading contribution given by 
l/m^, and the wave function renormalization with Zq ^ 1. By taking into account the 
dispersive effects, and the fact that ei± ■ e2± — 1 due to the small opening angle between 
ki and k2, one has 




{ki ■ k2){ei± ■ e2±) ~ - 

{ki ■ e2±)ik2 ■ £i±) ~ -uJiUJ2sm^9i2 ~ -3^1-^1 ujiuj2 ■ (72) 




In the case of weak magnetic fields B <^ Bcr and/or small energies u <^ m, one can see 
from Eq.(72) that the second term in parenthesis in Eq.(71) is sub- leading with respect to 
the first one and so it can be safely neglected. Finally, by using Eqs.(52), (61), (71), (72) 
and integrating over all the phase space, we find the following result for the absorption 
coefficient rfy-* induced by the massless mode 




rff --— ^ -- / diUi dU2iiUliU2f5iiU-iU^-iU2). (73) 



7 <-,^ , U , /,,,■* ,., /Q 
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Performing the last integral in the phase space and eliminating the delta-function, we get 

5 




-^^) , (74) 

38407r/ ' ^ ^ 



where Bcr = /i^ and /i = \/mK. The above result holds only provided that 

BcJ \m) 

Regarding the photon splitting in the case of a scalar coupling with an external electric 
field, the same results can be easily obtained from the above expression. At the leading 
order in the A/m^ expansion, the expression for the absorption coefficient is the same as 
in Eq.(74), provided that the external magnetic field is replaced with the electric one, and 
the allowed polarizations are 7|| -^ 7|| 7|| . 

Now we consider the case of a pseudoscalar coupling. If we neglect for the moment 
the effects of the small angles in the final momenta, it is not difficult to see that due to 
the parity violating coupling, the polarizations of the final photon states will be mainly 
opposite, in particular 7|| -^ 711 7_|_ or 7|| -^ 7_l /ii • ^'^ ^^is case, by taking into account 
the dispersion effects in the refraction indices and kinematics, we see that the condition in 
Eq.(62) cannot be satisfied, therefore the corresponding photon splitting is not allowed in 
this case. On the other hand, the process 7|| -^ 7|| 7|| , due to the pseudoscalar coupling, 
is forbidden. Same results hold for the case of an external electric field. 

It is interesting to compare the result of the absorption coefficient induced by scalar 
interactions in Eq.(74) with the corresponding one due to vacuum polarization effects in 
QED. In this last case one has [21] 

<^r-0.l(5^)°(^)%,n-. (76) 

where B^^^ = m1/e = 4.41 x lO^Tesla, with rrie and e the electron mass and charge 
respectively. If one considers as an example the values of scalar mass and coupling around 
m = lO^'^eV and A = lO^GeV, as for instance suggested by the central values of PVLAS 
data [17], fi = IMeV, which is incidentally quite close to the scale of the critical magnetic 



field in QED, namely yBcr . Then, by using Eq.(74), one gets 

4°) 10 ^rn ^ '' 



jQED 



From Eq.(77) follows that the QED contribution to the absorption coefficient is smaller 
than the scalar /pseudoscalar one, provided that /i < 5 x mg . We stress that this upper 
bound is stronger than the one obtained by requiring {n± — 1)'^^^ /{n± — 1) < 1, implying 
that /i < 24 X nie ■ However, for the value /x = IMeV, as suggested by PVLAS data, the 
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ratio in Eq.(77) is -qfen — 4 x 10^ which shows that the massless mode induces quite a large 
effect in the photon sphtting with respect to the QED one. It is clear that for laboratory 
magnetic fields of the order of 1 Tesla and B ^ B^r also the scalar/pseudoscalar contri- 
bution to the absorption coefficient induced by the massless mode is quite small, being 
suppressed by six powers of the external field. However, although strongly suppressed, 
this effect could have interesting application in astrophysical context, where for instance 
the intensity of magnetic field in the core of a supernova could be very large. 

The photon splitting through the massive mode remains the main mechanism for 
laboratory experiments, provided the critical condition ^ < 1 for external fields is satisfied. 
Indeed, by inspecting the ratio between the absorption coefficients induced by the massless 
(d^^) and massive mode (rf^), for external fields far above the critical value ^ <C 1, one 
has 



rff 




(78) 



Then, from Eq.(78) it follows that, ii fx> \/f^, the mechanism of photon splitting by 



2m' 

means of massive mode is dominant, and, for magnetic fields satisfying the condition 
B < SOrrig (^), it is also larger than the corresponding QED effect. 

Finally, we directly compare the QED background to the photon splitting process 
induced by the massive mode, when the magnetic fields are far above their critical value 
,^ ^ 1. In this case one gets, 

s!^ . 5.4 X 10-- (^y (^y (^) f^x) (By 

d^ \mj Km^J Km^J \m J VlTesla/ 

In conclusion, for laboratory magnetic fields of the order of Tesla, scalar/pseudoscalar 
masses and couplings in the range 10~^ eV < m < 10^ eV and 10^ GeV < A < 10^° GeV 
respectively, and photon beam energy in the range of 1 eV < uj < 10^ eV, the dominant 
mechanism of photon splitting is through the massive mode, and the QED background is 
negligible. 



5 Numerical results 

In this section we present the numerical results of a model independent analysis based on 
the photon splitting process in two photons mediated by scalar/pseudoscalar interactions. 
In particular, we will show that the search for photon splitting process by using laser 
experiments in the optical and X-ray frequency range, could allow to explore extensive 
regions in the {m — g) plane, depending on the intensity of the external magnetic field, 
where g = l/A^p. Remarkably, in the case of pseudoscalar (axion) coupling, these regions 
would be complementary to the ones already excluded by current laser experiments. 
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Figure 3: Upper bounds on m in tlie scalar case, versus tlie piioton energy to in eV, for three 
values of external magnetic field B = 0.1, 1, 10 Tesla, respectively from bottom to top. 

We will consider only the effect of the photon splitting induced by the massive mode, 
since, as shown in section 4.2, it is the dominant effect for laboratory experiments. We 
emphasize here that, fixing the values of the photon energy and external field, the critical 
conditions in Eqs.(37), (38) set strong upper bound on the scalar/pseudoscalar masses 
m. To get a feeling with numbers, in Fig.3 we show the results for the allowed values of 
scalar masses m versus the photon energy and for several values of magnetic field. We 
consider three representative cases of laboratory magnetic fields, namely |B| = 0.1, 1, 10 
Tesla. For example, for photon energies from optical up to the soft X-ray range, uj = 
(1 -200) eV, as for instance the free electron laser with very high peak brilliance from 
UV and (soft) X-ray sources at DESY and SLAC [18,19], and magnetic fields between 
1-10 Tesla, masses m up to 100 eV can be explored. However, when the photon energy is 
UJ ':^ m, y\B\, the mass upper bound is much below the corresponding to. 

For example, it is possible to realize in laboratory high energy photon beam by laser 
back-scattering on a primary electron beam [27], as in gamma-gamma colliders [28]. How- 
ever, notice that even with an energy beam of cu ~ 10^ GeV, which is beyond the capability 
of planned collider experiments, and a magnetic field of the order of 10 Tesla, one can 
explore massive modes up to m of the order of MeV. 

Now we will restrict our analysis to the case of low energy photons, in particular be- 
tween the optical and soft X-ray range, where the relevant scalar/pseudoscalar photon 
couplings can only be to photons and neutrinos. For neutrino masses of the order of meV 
the J ^>- up conversion process would be kinematically allowed and the F appearing in 
Eq.(51) should be identified with the total width. Then, the formula for the absorption 
coefficient of the 7 — > 77 process should be multiplied by the corresponding branching 



27 



lO"" 



10"- 





10 


<—• 




> 
O 


10 
10 


D) 


10 



10"' 



10"' 



CO = lev 




co=10eV 



co = 200eV ^ 



10"' 



m[eV] 





lO"-" 




10-* 




10"= 


<—• 




> 
O 


10"= 
10"^ 


O) 


10"= 




10-' 




10"' 



10"' 



co = 1eV 



(n = 10eV 




(0 = 200 eV 



10"' 



10" 10^ 

m[eV] 



10' 



Figure 4: Excluded areas at 95 % C.L. in the g = 1/A and m plane, corresponding to photon 
energies u; = 1 eV (blue), u; = 10 eV (green), uj = 200 eV (red). We have assumed a magnetic 
field of 5 Tesla and L = 10 m, a photon beam with momentum orthogonal to the direction of 
magnetic field {9 = vr/2) and number of photons per second dN^/dt = 10^^s~^. Results for an 
integrated running time of one day and one year are shown in left and right plots respectively. 

ratio. However, in the case of a PQ axion, the effects of axion couplings with neutri- 
nos on the total width can be neglected due to the fact that axion-neutrino coupling is 
suppressed by neutrino masses. Indeed, the interaction Lagrangian in this case would be 
L = gA^^lb^ ^A where 99^ is the axion field and qa oc m^/ K, while A should be identified 
with the coupling Ap appearing in Eq.(2). 

Since we would like to generalize the results of our analysis to the PQ axion case, 
we restrict ourselves to the case where scalar/pseudoscalar particles are only effectively 
coupled to two photons. The total width will coincide in this case with the width in two 
photons given by 



m 



647rA2 



(80) 



The same expression, as a function of mass m and coupling A, holds for both scalar 
and pseudoscalar interactions in Eqs.(l) and Eqs.(2). The fact the scalar/pseudoscalar 
particles are only coupled to photons greatly simplifies the analysis. In particular we can 
see that the expression for ^ does not depend on the coupling A and it is given by 



e 



1 



m 



12 



1024 7r2|k|4|B|4sin^^ 



(81) 



2 2 
~ a; — m . 



where |kp 

Now we consider a realistic experiment in which a monochromatic photon beam of 
frequency uj is traveling through a constant and homogeneous magnetic field of length 
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L. If we choose the direction of polarization of the magnetic (electric) field component 
of the photon beam to be parallel to the external magnetic field, 9 = 7r/2, we select 
the maximum coupling for scalar (pseudoscalar) contributions. The main signal consists 
in photon pairs with different energies Ui and uj2, restricted by the condition of energy 
conservation uj = uji + 002- In the case in which uj ^ m, the two photons are produced 
both almost forward and parallel to the original momentum of the photon beam, with a 
small opening angle 6 ~ O{m/uo) <^ 1.^ Then, excluded regions on the {m — g) plane 
can be set, for instance, by requiring that no significant number of events are observed 
at 95% confidence level. Since this process has practically no background, due to the 
fact that in QED this effect is very suppressed, the corresponding significant number of 
standard deviation associated to a number of observed events A^^ would be of the order 
\/Ns. In the particular case of 95% C.L. this implies Ns — 4. Then, the requirement 
that no significant number of events are observed at 95% C.L. would imply 

d; < -^ (82) 

where d^ is given in Eq.(51), with F and ^ given by Eqs.(80) and (81) respectively. For 
N^ we have assumed a representative laser brilliance corresponding to dN^/dt = 10^^/sec, 
as for instance in the case of free electron lasers. 

Excluded regions in the {m — g) plane corresponding to the upper limit in Eq.(82) 
are shown in Fig. 4 for a magnetic field of 5 Tesla and for different values of photon beam 
energy, namely cj = 1, 10, 200 eV. Here in the left (right) plot we show the exclusion 
regions corresponding to one day (year) of running time, respectively. As we can see 
from these results, by increasing the photon energy, larger values of A can be probed. 
In particular, with laser frequencies in the optical range 00 = 1 eV, values of A < 10^ 
GeV can be explored after one year of running, while A < 10^ GeV can be reached with 
UJ = 200 eV. The range of masses that can be explored with a 5 Tesla magnetic field, are 
just limited by the photon energy, namely m < u, for u = 1, 10 eV. However, for this 
value of the magnetic field and uj = 200 eV, the resonant region (^ ~ 1) is achieved below 
the kinematical upper bound m < u, and so for u = 200 eV only the range of masses 
up to m < 125 eV can be explored. Around the region close to the critical point, the d^ 
becomes insensitive to A as shown in the previous section, see Eq.(58). This is the reason 
why the shape of the lower part of the red area, corresponding to uj = 200 eV, near the 
end point of m is different from the other cases. The narrow width characterizing the end 
point region with large values of A, is of the order of F. 

^In detecting the signal, one could use a particular device where all the photons with energies very 
close to the primary photon beam energy to can be absorbed. Then, by placing a photon detector around 
the magnet, the split photons could be detected. It is beyond the scope of the present paper to analyze 
an efficient way for measuring the photon splitting, and this work should be considered as a theoretical 
proposal. 
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Figure 5: Results of right plot in Fig. 4 embedded in the g = 1/A and m plane with exclusion 
regions coining from laser experiments BFRT [15] and PVLAS [17]. Exclusion regions coming 
from solar experiments searches for axion from the sun (Solar CAST) [29] and the horizon- 
tal branch stars (HB stars) [30] are also shown. The inclined lines (green) correspond to the 
predictions in DFSZ [31] and KSVZ [32] axion models. See text for further details. 

Let us now discuss the consequences of these results for the axion physics. We remind 
here that the mass of the axion field connected to the PQ symmetry is related to the mass 
and decay constant of the pion by [4] 



niA 



z^l^lil + z) m^ U/JA = 0.6 meV x (lO^OGeV//^ 



(83) 



where Ja ^ 247 GeV is the scale where this symmetry is broken, /^ and m^ are the 
decay constant and mass of the pion and the current quark mass ratio z = mu/md. The 
couplings of the axion to SM particles are not only functions of the scale /a, but also 
model dependent. For instance, the constant /a is connected to the photon coupling 
gA'y = Ap^ in Eq.(2) by axial anomaly, and it is given by 



9A-y 
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31 + 2 



(84) 



where E/N is the ratio of electromagnetic over color anomalies and it is model dependent. 
The most popular models are DFSZ [31], where axion is embedded in Grand Unified 
Theories, and KSVZ axion is the so called "hadronic" class of axions [32]. In both these 
models the couplings gA'y are predicted in terms of Ja and are slightly different. The 
most stringent constraints on qa^ comes from cosmological and astrophysical arguments. 
However, there is a general class of experiments with laser that aim to produce axions in 
laboratory. They can be divided in two classes, one connected with photon-regeneration 
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and another one that analyzes the indirect axion effects on hght propagation. The first one 
is based on the Primakoff process [12], where a photon is converted in external magnetic 
field and then axions can be re-converted in photons with same mechanism after passing 
through a wall where all primary photons are stopped. The latter one probes the axion- 
coupling by measuring the change in the polarization of the photon after passing through 
magnetic field. 

In Fig. 5 we report the comparison between the present excluded areas in the {m — g) 
plane in the case of axion searches, by embedding the results contained in the right 
plot of Fig.4 corresponding to one year of integrated running time. The predictions of 
the models of DFSZ and KSVZ axions are also shown. The exclusion regions based 
on the solar experiments for axions from the Sun [29] (indicated with Solar (CAST) = 
CERN Axion Solar Telescope) and the constraints from horizontal branch (HB) stars 
[30] arising from a consideration of stellar energy losses through axion production are 
reported. Recently, the PVLAS experiment [17], reported the observation of a rotation 
of the polarization plane of light propagating through a transverse static magnetic field. 
This is shown as a point with error bars on the border of laser area in Fig.5. Notice 
that the signal seen in [17] is incompatible with the exclusion area from experiment [15]. 
Constraints coming from the galactic dark matter experiments [33] are not reported in 
Fig.5. However, this class of experiments strongly constrain only a narrow region in the 
{m — g) plane corresponding to axion masses in the range 10~^eV < rria < lO^^eV. We 
have also included the results of a recent proposal [18, 19] of laser experiment based on 
high-energy photon regeneration mechanisms, which is the X-ray analogous of the optical 
one. The free-electron lasers (FEL) at DESY's TESLA and dipole magnets of the type 
used in DESY's electron-proton collider HERA have been employed. The two horizontal 
lines correspond to the experimental analysis based on 1-1-1 and 200-1-200 HERA dipole 
magnets. 

As clear from the plot, the new areas in the (m — g) plane that can be probed by 
searching for photon splitting in two photons would allow to cover complementary areas 
not achievable by other two classes of laser experiments. It is worth noticing that both 
the predictions of KSVZ and DFSZ models fall inside the area that could be probed by 
the u = 200 eV laser experiment proposed here. As we can see from these results, the 
astrophysical constraints rule out a large region of the parameter space including the one 
we are interested in, as well as the region covered by laser experiments. However, recently 
it has been shown that in certain models these astrophysical constraints can be evaded 
[34]. 

Finally, we would like to emphasize that in the case in which a significant number of 
events of photon splitting 7 -^ 77 should be found, the values of the scalar /pseudoscalar 
mass and width can be easily reconstructed. This would allow to precisely tune the 
external magnetic field and photon energy in order to approach the resonant region, 
where the effect of photon splitting could be largely amplified. Regarding the origin of 
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the interaction, if there is a scalar or pseudoscalar couphng, this can be disentangled by 
analyzing the polarization of the final photons with respect to the direction of the magnetic 
field or by analyzing their angular distributions, provided that a sufficient number of events 
are observed. 

Let us also shortly discuss the case of the Higgs boson effects in this context. We 
start with the case of high energy photon beam traveling in a constant magnetic field and 
with the coupling to a light Higgs boson of mass rriH — 120 GeV. The main decay mode 
of the Higgs in this scenario is in bb. In the effective propagator there would be a very 
light mode, and a massive mode of the size of the Higgs mass. In order to generate the 
'-f ^ bb decay in external field one needs first to check if the condition (36) is satisfied 
for the corresponding massive mode. As we discussed in the previous sections this would 
require too high critical magnetic fields, when the energy of the photon beam is of the 
order of TeV, which cannot be realized in the laboratory. Nevertheless, there might be 
the possibility of a sizeable effect on the photon splitting mediated by the Higgs effective 
coupling in the framework of astrophysics, for example in the core of a supernova, where 
magnetic fields are very large. However, in that case the approximation of considering 
only constant and homogeneous magnetic fields is not correct, since also plasma effects 
should be taken into account. When magnetic fields are inhomogeneous, momentum can 
be absorbed and the Higgs boson could in principle be produced on-shell by means of 
the Primakoff process. Although the analysis of this latter effect in astrophysical context 
should be quite interesting, it goes beyond the scope of the present paper. 

6 Conclusions 

In this paper we have investigated the analytical properties of the effective photon prop- 
agator in an external homogeneous and static magnetic or electric fields, in the presence 
of scalar /pseudoscalar couplings. These results have been obtained by summing up in the 
photon propagator the relevant class of Feynman diagrams. They include the ones with 
mixing term of photon with scalar/pseudoscalar fields, where in the scalar /pseudoscalar 
propagator the corresponding self-energies have been exactly summed up. Then, we an- 
alyzed the solutions of the associated dispersion relations. 

Due to the presence of an external field, the standard photon dispersion relations in 
vacuum will be modified. While the effect of the real part of scalar/pseudoscalar self- 
energy can be re-absorbed in the corresponding mass term of the scalar/pseudoscalar, the 
imaginary part, connected to the corresponding width F, could play a crucial role in the 
photon dispersion relations giving rise to a non-vanishing contribution to the imaginary 
part of the effective photon propagator. In other words, the presence of F can induce a 
non- vanishing contribution for photon absorption coefficient. 

As known, two new propagation modes are allowed for the photon, in addition to 
the usual one which is not affected by the external field contributions. The lightest one 
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is associated to a light-like mode with refraction index n > 1, while the heaviest one 
corresponds to a massive mode with mass of the order of the scalar /pseudoscalar one. 

The new aspect of our work with respect to previous studies concerns the inclusion of 
the effects of the imaginary part of scalar /pseudoscalar self-energy in the effective photon 
propagator. In particular, we have shown that the massive mode can be allowed only when 
the contributions of the scalar /pseudoscalar width F are smaller or comparable to the 
mixing effects, see Eq.(36). This will give rise to a critical condition for the external field, 
which is absent in the limit F — > 0, depending on the photon energy, scalar /pseudoscalar 
masses and the two-photon coupling. A finite width for the photon can then be induced 
for the massive mode, being proportional to F, generating a non-vanishing value for the 
photon absorption coefficient d^. 

Although very small, being suppressed by the scalar/pseudoscalar width, the d.y can 
be sizeably enhanced due to a strong resonant phenomenon. In particular, we have found 
that a potentially large contribution to the splitting conversion 7 ^ 77 could be generated 
from the scalar/pseudoscalar decay in two photons when the external field approaches its 
critical value. 

We have also analyzed the contribution to the absorption coefficient induced by the 
massless photon mode. In the case of an external magnetic field and a scalar interaction, 
we find that the probability of photon splitting in two photons turn out to be very small, 
being suppressed by (B/Bcr)^ where (in relativistic unities) Bcr = mA. Nevertheless, this 
effect, depending on mass m and coupling A, could be much larger than the one induced by 
the QED vacuum polarization in the presence of an external magnetic field. For instance, 
for characteristic values oi m = 10~^eV and A = lO^GeV, we find that d^ is about five 
order of magnitude larger than the corresponding QED one. The leading contribution to 
(i-y is provided by the small effects of dispersions in the refraction index. For the case of a 
pseudoscalar interaction, due to the parity violating coupling and kinematic factors, the 
photon splitting induced by the massless mode is forbidden. 

We have analyzed the consequences of these results for a new kind of laser laboratory 
experiments, based on the searching of the photon splitting 7 ^ 77 conversion in external 
constant and homogeneous magnetic field. In particular, we have considered the photon 
splitting phenomenon induced by the massive mode. By taking the case of high brilliance 
lasers with dN^/dt = 10^^/sec, in the range between optical (cj ~ 1 eV) and the low 
X-ray {00 ~ 200 eV) frequencies, magnetic fields of 10 meters long and of the order of 
5 Tesla, we show that it is possible to probe a large area of the m and g = 1/A plane, 
provided that the process 7 — > 77 can be efficiently detected. The probed areas for 
scalar and pseudoscalar case are slightly different. We find that in the case of a; ~ 200 
eV, scalar/pseudoscalar masses can be probed up to m < 100 eV , while in the case of 
uj = 1 eV the region up to the maximum value allowed by the energy conservation, i.e. 
m < 1 eV can be explored. Moreover, the sensitivity on the scale A associated to the 
two-photon coupling can become , in one year of running time, close to A ^ 10^ GeV and 
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A ~ 10^ GeV, corresponding to the case of cj = 1 eV and u = 200 eV laser frequencies, 
respectively. 

We have also compared our results with the present bounds from axion searches. In 
particular, by restricting our predictions to the case of pseudoscalar couplings, we found 
that a large area in the m and g plane of the axion can be tested. Remarkably, this 
area is complementary to the ones already explored by the present laser experiments, 
like PVLAS and BFRT, and to the sensitivity area of new recent proposal based on the 
photon-regeneration class of experiments with X-ray free electron laser facility at HERA. 
We have also shown that some predictions of the KSVZ and DFSZ models could be tested 
by means of 7 — > 77 searches in external magnetic fields. 
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Appendix A 

In this appendix we report the exact solutions of the mass-gap equation (26). This is a 
cubic algebraic equation and can be analytically solved by using the known cubic formula. 
In particular, by substituting k'^ = m? + e in Eq. (26), this can be simplified as 

x^ + x' + x (er - a) + ^r = , (85) 

where x = e/m^, C,r = ^'^/iri^ and ^^ = A/m^. Notice that in writing the gap equation we 
have identified the ImS(A;^) with —mT. In general this statement is not correct since the 
imaginary part of scalar/pseudoscalar width depends on /c^, in particular for the photon 
splitting it is proportional to (A;^)^/A^, see Appendix B for more details. However, the 
approximation to set ImS(A;^) = —mT is valid only when /c^ ~ 0{m'^), that is for the 
massive solutions. On the other hand, for the massless mode where, /c^ ~ 0{A/m^), the 
ImS(A;^) can be neglected with respect to the m'^ term and one can safely switch off the 
ImE(/c^) in the scalar/pseudoscalar propagator provided that A/m'^ <^ 1. Therefore, in 
order to simplify the problem, one can solve the gap equation by using ImS(A;^) = mT 
and then, in order to recover the correct result for the massless mode, set ^r -^ on the 
corresponding massless solution. 
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By using this approach, we solve the cubic equation above and find the following 
expressions for the solutions of Eq. (26), namely k"^ = Mq and k"^ = M]_ with 

M^ = lim m^ f I + X+ + X_ 
" 5r--o V3 

Ml = m^(---{X^ + X^)±t—{X+-X^)], (86) 



where 



and in our case 



X± = (i? ± Vof^^^ , D = Q^ + R'^ (87) 



g = ^ , i?= ^^ . (88) 

The massive solutions k'^ = M|. are real only when the condition Z^ < is satisfied, where 

D = ^{-e^ + A^r-^eA + 20U^r + l2 (^i^r - er^) + 8^^ + 4^?} . (89) 

In our problem, the parameters ^a <^ 1 and ^r <^ 1 are very small quantities. We 
remind here that these are proportional to inverse powers of the effective coupling A with 
two photons, which is assumed to be much larger than the typical scales appearing in 
the numerator of ^r and .^Aj namely the photon energy u, momentum k, and external 
field. From a first inspection of Eq.(89), one can see that, in order to have real solutions 
(-D < 0), these parameters should satisfy the following hierarchy ^r ~ Ca- This is also 
consistent with the fact that .^r and ^a scale as ^r ~ 1/A^ and .^a ~ 1/A^ respectively 
as a function of A. In particular, by dropping the terms of order 0{C,%) in Eq.(89), the 
condition D < is equivalent to 

^A > a^°^ = 2^ (90) 

which is exactly the inequality in (36) derived under the approximation of neglecting 
terms of order ^^. 

There is only one solution of D = 0, corresponding to real and positive values of ^a 
and ^-p, and it is given by 

f - f 2 2(1- 27eA) , HjU) .... 

where the expression for the function H{x) is 

H{x) = (8 + 540x - 729a;2 + 3 v^S v^ (8 + 27a;)^/^) . (92) 
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As we can see from these results, there always exists a critical value for the external fields 
for which the massive solutions become degenerate, implying the appearance of a double 
pole in the effective photon propagator. 

Now we will provide below the corresponding solutions obtained by expanding the 
exact ones in terms of the small parameters .^a and .^r- Let us first define ^^ being the 
value of .^A satisfying the exact solution of Eq.(91). Then, by expanding ^a in powers of 
6 = r/m, we obtain, up to terms 0{6^), the following expression 

Ca = Ca^'^ (l + Ri5) ) (93) 



where 

R{x) =x-lx^ + x^ -^x'^ + 8x^ + 0{x^). (94) 

Finally, the corresponding next-to-leading order corrections, in ^a expansion, to the 
solutions in Eqs.(30) are given by 

Mo = -^(l-U + ^^A 

Ml ^ ,„^{l4(l±;i-5j-^^l±5^j| (95) 

where the symbol ^ = 4^r/'CA i^ ^^e same as appearing in Eq.(31). Notice that for the 
solutions in (95) the range < .^ < 1 is not valid anymore, and, according to Eq.(93), 
must be modified as < .^ < l — 6 + 0{6'^). It can be easily checked that for the maximum 
value ^"^'^^ = I ~ S, the difference M^ — M^ vanishes up to terms of order ^(^a)) within 
the validity of the approximated solutions. 

Regarding the massless mode, the expression for Mq can be easily re-obtained by 
setting to zero the scalar/pseudoscalar width F directly in the gap equation. This will 
reduce the cubic equation in Eq.(26) to a second order one. In particular, by setting 
ImE(/c^) to zero in Eq.(26), one gets a more compact solution for the massless mode, 
equivalent to the first solution in Eq.(86), given by 



m 



M', = '-^i^l-^l+4U) , (96) 

where the leading orders in ^a expansion easily recover the corresponding results in 
Eqs.(30) and (95). 

Appendix B 

In this appendix we provide the calculation at 1-loop of the imaginary part of scalar 
(pseudoscalar) self-energy T,s{p), induced by the interactions in Eqs.(l),(2). While the 
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real part of T,s{p) is divergent, its imaginary part is finite at 1-loop and can be easily 
calculated by making use of the optical theorem. In particular, the Ira[T,s(p){k'^)] is 
connected to the square modulus of the off-shell scalar/pseudoscalar decay in two on-shell 
photons, which schematically can be represented as 

1 
2 



Z Im 



.. ■ . . _____^ 



<:> 






(97) 

where Z is the renormalization constant of scalar/pseudoscalar field, the continuous and 
dashed line stands for the scalar/pseudoscalar and photon fields respectively, and k, ki^2 
are the corresponding momenta. In the left-hand side of (97), the vertical dashed lines 
correspond to the Cutkosky rule which selects the imaginary part of the photon propaga- 
tors in the loop, by replacing in each of them l/(p^ + ie) -^ —7t6{p'^)6{uj), with p^ = PfiP'^ 
and uj = po, where p^ indicate the generic 4-momenta of photon lines in the loop. The 
factor 1/2 in front of the amplitude square is due to symmetry factor connected to the two 
photon lines. By using the Cutkosky rule and the fact that at this order in perturbation 
theory Z = 1, one obtains 

where V^^^^^\ki, k2) is the two photon vertex for the scalar (pseudoscalar) field given in 
Eq.(55), and the ki^2 are the corresponding final photon momenta on shell, /c^g — 0- 
The sum X^poi is extended to all photon polarizations. Notice that, as an approximation, 
we have not considered the effect of the external field corrections, and so the photon 
dispersion relations are not modified. 

Finally, by using the results in Eq.(55), we get, for the contribution at 1-loop to the 
self-energy. 

When Eq.(98) is evaluated on-shell the scalar (pseudoscalar) width ^s{p) on-shell is recov- 
ered by means of the relation Im[S(A;^)] 1^2=^2 = —mT, in particular 

in 



which is in agreement with Eq.(80). 
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